Charge-density- wave formation by Van Hove nesting in the a-phase of Sn/Ge(lll) 



J. Gonzalez 

Instituto de Estructura de la Materia. Consejo Superior de Investigaciones Cientificas. Serrano 123, 28006 Madrid. Spain. 

(February 1, 2008) 



O 
O 
O 

(N 
< 



X3 
O 

o 



> 

o 

o 
o 
o 



o 



X 
S3 



We study the role of electron correlations in the formation of the surface charge-density-wave state 
in the Sn/Ge(lll) interface. The Fermi energy of the overlayer is treated as a dynamical variable, 
which undergoes a substantial renormalization by the interaction. We show that the Fermi level 
turns out to be pinned to a Van Hove singularity in the density of states, which explains the formation 
of the charge-density-wave, the observation of a very flat band in photoemission experiments and 
the reduction of the spectral weight in the low-temperature phase. 



The formation of states with rearrangement of the elec- 
tronic charge is accompanied in general by interesting 
properties in condensed matter systems. A recent ex- 
ample has been provided by the study of interfaces like 
Pb/Ge(lll) §|, Sn/Ge(lll) |i or K/Si(lll) §. All 
of them have in common the existence of a phase transi- 
tion with the formation of a surface charge-density-wave 
(CDW) in the low-temperature phase. In the first two 
cases the adsorbate phases have at room temperature the 
same (-\/3 x -\/3)i?30° structure, and they both show a re- 
duction of the symmetry to a (3 x 3) ground state below 
the critical temperature. However, while the Pb/Ge(lll) 
overlayer becomes insulating in the charge-ordered state, 
the Sn/Ge(lll) remains metallic below the transition. 
In the first case, the analysis of the surface band has 
revealed nesting along significant portions of the two- 
dimensional Fermi surface which would account for 
the formation of the CDW. In the a-phase of Sn/Ge(lll), 
however, no similar effect seems to be at work ||,^, and 
the reason why the ground state of the system is led to 
reorder its charge is not fully understood yet. 

It is believed that correlation effects may play an im- 
portant role in the development of the charge-ordered 
state [||j7|. The relatively small surface bandwidth of 
these systems (~ 0.5 eV) leads to consider the on-site 
Coulomb repulsion as one of the relevant interactions. It 
has been also suggested that a nearest-neighbor repulsion 
may play an important role in producing the 3x3 surface 
periodicity 

In this letter we develop the idea that the electron cor- 
relations are essential to understand the low-temperature 
phase of Sn/Ge(lll), while relying on the fact that no 
strong-coupling effects need to be invoked for that pur- 
pose. This is plausible as long as the overlayer remains 
metallic below the transition. We will show that the 
(\/3 X \/3) phase is unstable towards the formation of 
saddle points near the Fermi surface in the electronic 
dispersion relation. This is in agreement with the obser- 
vation of a very flat band below the Fermi level in the 
angle-resolved photoemission experiments carried out in 
the CDW phase M. The Van Hove singularity controls 



the dynamics of the system since, as we will see, this finds 
energetically favorable to readjust the bands in order to 
place the Fermi surface closer to the saddle points. As 
well as nesting of the saddle points near the Fermi surface 
gives rise to a spin-density-wave instability in the repul- 
sive t — t' Hubbard model [p|-p^ , it is able to explain the 
formation of the CDW under the assumption of a more 
general interaction in the present instance. 

We will take an extended Hubbard model with on-site 
U and nearest-neighbor V repulsive interactions on the 
triangular lattice of the Sn/Ge(lll) overlayer to illus- 
trate these ideas. The adatom dangling bonds give rise 
to a band that reaches a minimum at the corner of the 
hexagonal (\/3 x ^/i) Brillouin zone, and a maximum 
close to the center . We use a tight-binding approxima- 
tion with nearest neighbor hopping t and next-to-nearest 
neighbor hopping t' to parametrize the Fermi line of the 
(V3 X -\/3) phase as given in Ref. §. Although the pre- 
cision of the fit to the dispersion relation is limited with 
this approximation, we are able to capture in this way the 
main topological features, which play an essential role in 
the following considerations. 

The topology of the dispersion relation is actually un- 
stable under the smallest periodic lattice distortion. This 
has to be taken into account, as it has been recently ob- 
served that one of the Sn atoms in the (3 x 3) unit cell 
is slightly displaced upwards with respect to the other 
two [ |l3| , |T^ . By taking a hopping parameter equal to 0.9t 
between inequivalent Sn atoms in the (3 x 3) unit cell, 
we already observe the appearance of two subbands of 
very different character. The original valence band for 
the (\/3 X a/3) lattice splits along the FK direction of the 
(3x3) Brillouin zone. Thus, there are an upper and a 
lower subband that join at the apex of the conical dis- 
persion shown in Fig. |l|, characteristic of the underlying 
honeycomb lattice, and a middle subband that develops 
a saddle point at the boundary of the Brillouin zone, as 
shown in Fig. ||. The density of states redistributes itself, 
having a linear dependence in energy near K in the up- 
per and lower honeycomb-type subbands while becoming 
divergent at the M point in the middle subband. 
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FIG. 1. Lower subband of the (3 x 3) phase. The upper 
subband (not shown) displays symmetric cusps at the six cor- 
ners of the (3 X 3) BriUouin zone. 




FIG. 2. Middle subband of the (3 x 3) phase. The saddle 
points lie at the boundary of the (3 x 3) BriUouin zone. 

Before switching the interaction, the Fermi level has a 
nominal location below the mentioned Van Hove singu- 
larity. However, the relative position of the Fermi level 
in the different subbands of the model may suffer signifi- 
cant changes by the effect of the electron correlations near 
the singularity. The forward scattering vertex Ffwd giv- 
ing the four-point interaction at small momentum trans- 
fer k gets corrected by the iteration of the particle-hole 
polarizability. For the saddle-point dispersion relation 
e(k) fa k^— ky, this object diverges in the static hmit as 
~ log(A/max{|e(k)|, |e_F|}), where ep is the distance of 
the Fermi level to the Van Hove singularity and A is a 
cutoff of the order of the bandwidth ||l^. In the case of 
a spin-independent interaction, this leads to a reduction 
of the vertex at zero momentum-transfer 
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where the bare couphng is iJfwd = '^(Pi, P2; Pi, P2) — 
f^(Pi7 P2; P2, Pi), in terms of the original four-fermion 
interaction, and c is the prefactor of the particle-hole po- 
larizability . This renormalization of Ff^d has impor- 
tant consequences, since it encodes the sum of dominant 
logarithmic contributions that give rise to the constant 
shift of the chemical potential in the self-energy correc- 
tions In this respect, it has been shown that the 
Van Hove singularity has a natural tendency to pin the 
Fermi level in a model with bare repulsive interactions 
p8| , pj[ . In the present case, this can be understood in 
terms of the uneven screening of the forward scattering 
in the different subbands of the model, as we show in 
what follows. 

Let us model the problem of the interaction between 
the two subbands in Figs. |l|and ||by taking the divergent 
density of states near the saddle-point structure 
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and a vanishing density of states for the lower subband 
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The physical value of the Fermi energy comes from the 
balance between the nominal value n of the chemical po- 
tential and the shift of the one-particle levels to higher 
energies due to the repulsive interaction. The filling level 
of each subband is determined by the self-energy correc- 
tions, bearing in mind that the upward displacement of 
the levels depends in turn on the charge present in all the 
subbands. The filling level epi in the middle subband is 
given by the equation 



e_Fi = M - r 
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(4) 

where we have introduced a coupling constant gfwd that 
parametrizes the Coulomb repulsion exerted on the mid- 
dle subband by the charge present in the lower subband. 
Analogously, the filling level e f2 in the latter is computed 
in the form 



eF2 = A* - C/f 



wd 



(5) 

We remark that epi and eF2 are measured in the refer- 
ence frames in which the dependences rt*-^^ (e) and n*-^^ (e) 
are as given in Eqs. (^) and (^). Thus, the fact that epi 
and £i?2 are nominally different after renormalization is 
just an artifact of that convention. The physical picture 
is however the opposite, namely that the one-particle lev- 
els are shifted to higher energy by a different amount in 
each of the subbands, up to a point in which the respec- 
tive Fermi levels reach the common chemical potential. 
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The coupled set of equations (Q) and (||) gives rise to 
nontrivial physical effects, as a consequence of the non- 
linearities introduced by the divergent density of states 
n^^\e) and the renormahzation of FfwdCe) close to the 
Van Hove singularity. It is convenient to solve for the 
location of epi and Sf2 in terms of the total charge N in 
the system, given by 



iV 



(6) 



The most remarkable effect is that there is not a one-to- 
one correspondence between N and the respective filling 
levels epi and £^2 for the two subbands. The different 
branches of the solution are represented in Fig. || for 
the particular values f/fwd — 4. OA and ^fwd — 3. OA. The 
parameter [3 has been chosen equal to 3.0, and a has 
been set so that there is the same number of states in 
the lower and in the middle subband. 
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FIG. 3. Self-consistent solutions for the filling levels in the 
middle and lower subbands vs total charge in the overlayer. 

At low values of N , the filling of the middle subband 
with the Van Hove singularity proceeds in a regular way, 
with a monotonous increase of epi- There is a point, 
however, in which two other locations of epi become pos- 
sible, closer to the singularity in the density of states. In 
these instances, the corresponding filling level eF2 in the 
second subband suffers an abrupt decrease with respect 
to the expected value. It is interesting to discern what 
of the possible solutions is most favorable energetically. 
We have plotted in Fig. ^ the value of the total energy E 
versus the total charge N . We see that above a certain 
value, that is A^ ~ 1.4 in our case, any of the two fill- 
ing levels close to the Van Hove singularity give the most 
stable configuration of the system. This is in agreement 
with previous analyses of the pinning of the Fermi level 
of electrons near a Van Hove singularity ||l8|,|l9| . 
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FIG. 4. Total energy of the difi'erent solutions shown in 
Fig. 1^ vs total charge N (the symbols for a given solution 
match in the two figures). 

The result that turns out to be valid under very gen- 
eral conditions is the existence of an intermediate range 
of filling levels that are forbidden below the Van Hove 
singularity. For the corresponding values of the charge 
N ^ this finds more favorable to fill the Fermi sea up to the 
Van Hove singularity, at the expense of the charge in the 
other subband. In general, there is a critical value of N in 
which the filling level sfi jumps discontinously from the 
regular evolution upon doping to a position much closer 
to the Van Hove singularity. 

Coming back to the physics of the Sn/Ge(lll) over- 
layer, a description consistent with the plots in Figs. |l| 
and H leads to a nominal value of e fi that is about 25% 
off the Van Hove singularity, relative to the energy dif- 
ference between the level of the latter and the bottom of 
the subband. This nominal value falls certainly in the 
forbidden region, what leads us to conclude according to 
the above arguments that the Fermi level is pinned close 
to the Van Hove singularity developed in the system. 

This abrupt renormahzation of the Fermi level to- 
wards the Van Hove singularity is able to account for 
all the prominent experimental features found in the 
Sn/Ge(lll) interface, as we show in what follows. 

1. Strong correlations in the low-temperature phase. 

The different photoemission experiments have in com- 
mon the observation in the low-temperature phase of a 
lower band that disperses along the FM direction of the 
(3 X 3) Brillouin zone, together with another band with 
higher energy and much less dispersion [|jl3). The point 
is that this intermediate band is found at an energy sensi- 
bly lower than predicted by LDA calculations . Quite 
remarkably, all the points along the FM direction are 
found below the Fermi level. After taking into account 
the experimental overestimation of the binding energies 
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remarked in Ref . Q , the M point turns out to be close but 
below the Fermi level, contrary to the estimates carried 
out neglecting many-body effects. This is actually the 
experimental signature of the pinning to the Van Hove 
singularity we have discussed above. 

It is also clear that the width of the middle subband 
is sensibly smaller than predicted by LDA calculations 
P,p^ . This reduction in the bandwidth is also consistent 
with the renormalization of the kinetic energy of electrons 
close to a saddle point, which leads to a flatter shape of 



the dispersion relation |15 



2. Formation of the surface CDW. 

The saddle points of the middle subband are at the M 
points of the (3 x 3) Brillouin zone, so that they are sep- 
arated by momenta that correspond actually to the 
wavevectors of the CDW structure observed experimen- 
tally. This Van Hove nesting by momenta connecting 
the saddle points near the Fermi level leads to the CDW 
instability. 

In the context of an extended Hubbard model and with 
the inclusion of the electron-phonon interaction, it has 
been shown that it is natural to have a negative cou- 
pling f/Ki for the four-fermion interaction with momen- 
tum transfer Q (its bare value is actually U — QV, 
in terms of the parameters of the extended Hubbard 
model). We remark that, with a divergent density of 
states close to the Fermi level, the slightest attractive in- 
teraction {UKi < 0) triggers the CDW instability. The 
particle-hole susceptibility with momentum is itself 
logarithmically divergent, and the iteration of particle- 
hole diagrams leads to an effective vertex at momentum 
transfer K, 



rK.(£) 



l + {c'/{27r^A))UK, log|A/e| 
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Thus, the attractive interaction is overscreened and 
renormalized towards more intense attraction at low en- 
ergies. 

Under these conditions, the only response function 
that shows a divergence at some critical frequency corre- 
sponds to a charge modulation with wavevectors K.;. The 
gap that opens up at the saddle points can be actually 
estimated from the frequency at which the denominator 
in Eq. vanishes. The order of magnitude thus ob- 
tained (tens of meV) is consistent with the temperature 
of the transition to the CDW. 

3. Metallic behavior and depletion of the density of states 
in the low-temperature phase. 

The metallic behavior of the (3 x 3) phase follows nat- 
urally from our description of the different subbands in 
the system. Nesting of the saddle points at the Van Hove 
filling and the formation of the surface CDW lead to the 
opening of a gap and to a partial destruction of the Fermi 
line around the saddle points . The divergent density 



of states induces by itself a significant reduction of the 
quasiparticle weight, as appreciated in the computation 
of the self-energy corrections for the saddle point disper- 
sion relation This is consistent with the decrease 
of the photoemission intensity about the Fermi level ob- 
served in the CDW phase |4|. The lower band is less 
sensitive to these renormalization effects, and the metal- 
lic properties come mainly from the band crossing the 
Fermi level below the cusps shown in Fig. |^. 

In conclusion, the pinning of the Fermi level to a 
Van Hove singularity in the density of states leads to a 
phenomenology for the Sn/Ge(lll) overlayer consistent 
with the experimental observations of CDW formation 
and photoemission measurements, which are challenging 
more conventional theoretical approaches. Further pre- 
dictions of our approach should also be tested experimen- 
tally in Sn/Ge(lll), in particular the renormalization of 
the lower band and the geometry of the Fermi line around 
the corners of the (3 x 3) Brillouin zone. 
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